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1 Introduction 



It has been widely recognized that one-dimensional antiferromagnets with half- 
integer and integer spins have dramatically different excitation spectra. The 
original theoretical prediction by HaldaneJI| that Heisenberg chains with half- 
integer spin are gapless, whereas those with integer spin are gapped, has been 
confirmed experimentally ||. To gain insight into the physics underlying this 
beautiful result, one may study systems intermediate between spin S = 1/2 and 
S = 1. The simplest of these is the "Heisenberg spin ladder", which has isotropic 
couplings J|| along the chains and J± between them. Although such systems with 
both antiferromagnetic 0,0,0, @, an d ferromagnetic |8[l,H, [110|1 ,pT| inter- 
chain couplings have been the subject of considerable recent theoretical interest, 
certain problems remain unresolved leaving room for our contribution. 

In this paper we present our analysis of a weakly coupled spin ladder J\\ » 
| Jj_ | . We have found this limit more interesting from the theoretical point of view 
for it allows us to study the crossover between the gapless spin S = 1/2 regime 
and the strong coupling limit corresponding to the S = 1 chain. Unfortunately, 
our results have only qualitative validity for the known experimental realizations 
of double chain ladders ( Sr n _!Cu n+ i0 2n and (VO) 2 P20 7 []I3]]) where both 
exchange integrals are of the same order. 

The paper is organized as follows. In Section 2 we derive the continuous 
version of the spin ladder Hamiltonian for the case of identical chains. To achieve 
this we employ the bosonization approach, but the resulting effective theory is 
most simply represented in terms of fermions. In this representation the effective 
Hamiltonian of the spin ladder contains four species of weakly interacting real 
fermions^. Three of these modes comprise a degenerate triplet and the remaining 
one lies above having a mass approximately three times as big. The magnitude 
of the mass gaps is of the order of the interchain exchange. For any sign of 
the interchain coupling, the leading asymptotics of the correlation functions are 
determined by the triplet of Majorana fermions as for the S = 1 chain [fLl j. The 
fact that the low energy sector of the model is essentially a free theory makes it 
possible to obtain non-perturbative expressions for asymptotics of all correlation 
functions. This is done in Section 3. In Section 4 we discuss a situation where the 



1 The difference between ordinary (Dirac) and real (Majorana) fermions is that the latter 
ones have only positive energies e{p) = \/p 2 + m 2 . Therefore one can always describe one Dirac 
fermion as a superposition of two Majorana fermions. 
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ladder consists of inequivalent chains. It is shown that, in the limit when exchange 
integrals on the chains strongly differ, the low-lying excitations are described by 
the 0(3)-nonlinear sigma model. The adequacy of this treatment is guaranteed 
by the fact that this sigma model has a small bare coupling constant. As we have 
pointed out, the excitation spectrum of the spin ladder always has a gap. The 
appearence of such a gap is usually related to some symmetry breaking process. 
In Section 5 we discuss this process and derive the corresponding order parameter 
which happens to be nonlocal (the so-called string order parameter). The paper 
has a Conclusion and two Appendices where we provide technical details about 
bosonization and string order parameters. 



2 Coupling of identical chains; the Abelian 
bosonization. 

In this Section we apply the Abelian bosonization method to the spin-ladder 
model 

H = J \\ E E^'W-sA + ^ + ^ESiW-SsM (i) 

5=1,2 n n 

describing two antiferromagnetic (Jn > 0) spin-1/2 Heisenberg chains with a weak 
interchain coupling (| J±\ <C Jy) of arbitrary sign. Abelian bosonization is a well- 
known procedure, but for the sake of completeness we briefly overview it in the 
Appendix A. In the continuum limit, the critical properties of isolated S = 1/2 
Heisenberg chains are described in terms of massless Bose fields <j)j(x) (j = 1,2): 

#o = ^ E fdx[ nj(rr) + (d^x)) 2 ] (2) 

Z j=l,2 J 

where the velocity v s ~ J\\clq and IT,- are the momenta conjugate to (j>j. The 
interchain coupling 

H± = J±a J dx [ J\(x) ■ J 2(2) + ni(x) • n 2 (x) ] . (3) 

is expressed in terms of the operators 3j(x) and Hj(x) which represent, respec- 
tively, the slowly varying and staggered parts of local spin density operator and 



are defined in the Appendix A. According to (|106|) , the current-current term in 
(|3p is marginal, while interaction of the staggered parts of the spin densities is 
strongly relevant. So we start our analysis by dropping the former term (its role 
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will be discussed later). Using then bosonization formulas ( |105| ) for Tcij(x), we get 

H± = — / dx [- COS V / 27r(01 + 02) + 77 COS V / 27r(0l - 02 ) 

7T z ao J A I 

+ 003^(^1-^2)] 

where 3 -(sc) is the field dual to <f>j(x). Denote 

m = — — (4) 

Z7T 

and introduce linear combinations of the fields 0i and 02: 

01 ± 02 



V2 



(5) 



The total (0+) and relative (0_) degrees of freedom decouple, and the Hamilto- 
nian of two identical Heisenberg chains transforms to a sum of two independent 
contributions: 

H = H + + H_ (6) 
H+{x) = ^(nl + (d x <P + ) 2 )+—cosV^<P + (7) 

HJx) = - (n 2 _ + (9 X 0_) 2 ) + — 008^0- + — cos V^e^ (8) 

2 v ' 7ra 7ra 

In the above derivation, the Ji ■ J2-term has been omitted as being only 
marginal, as opposed to the retained, relevant rii ■ ri2-term. It is worth men- 
tioning that there are modifications of the original two-chain lattice model for 
which the Ji • J2-term does not appear at all in the continuum limit, and map- 
ping onto the model (H) becomes exact. In two such modifications, the interchain 
coupling is changed to 

^ ) = ^ESiW-[s 2 H-S2(n + i)] (9) 

or 

H[ B) = x DSi (n) - Sx(n + 1)] • [S 2 (n) - S 2 (n + 1)] (10) 

The structure of these models explains why the low-energy physics of two weakly 
coupled Heisenberg chains must not be sensitive to the sign of the interchain 
coupling J±. This conclusion is in agreement with recent results of Ref.7. 

Let us turn back to Eqs. (]?]) and @. One immediately realizes that the critical 
dimension of all the cosine-terms in Eqs. (0), (H) is 1; therefore the model (ra) is a 
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theory of free massive fermions. The Hamiltonian H + describes the sine-Gordon 
model at j3 2 = 4tt; so it is equivalent to a free massive Thirring model. Let us 
introduce a spinless Dirac fermion related to the scalar field <p + via identification 



Using 



iPr, l (x) ~ (27ra )- 1/2 exp (±iV^ <f>+- R , L ^)) ( n ) 
— cos VAtt <j>+(x) = i [ ip R (x)ipL{x) — h.c. } 



we get 



H+(x) = -iv s (ip R d x i/j R - i)\d x i) L ) + im(ip R ip L - ^\i\) R ) (12) 
For future purposes, we introduce two real (Majorana) fermion fields 

& = -^> ^ = —^~, {u = R,L) (13) 
to represent H + as a model of two degenerate massive Majorana fermions 

H+ = H m [Z\+HrM (14) 

where 

#m[£] = - y d x U - a + im Ml (15) 

(flm^J has the identical form). 

Now we shall demonstrate that the Hamiltonian H_ in (||) reduces to the 
Hamiltonian of two different Majorana fields. As before, we first introduce a 
spinless Dirac fermion 

X r,l{x) ^ (2na )- 1/2 exp (iiv 7 ^ <p- AL {x)) (16) 



cos V4n cf)^(x) = i [x r (x)xl(x) - h.c], 



ira 

— cos V4tt B-(x) = -i [ Xr(x)xl( x ) ~ h - c 



TTCto 

Apart from the usual mass bilinear term ("CDW" pairing), the Hamiltonian i7_ 
also contains a "Cooper-pairing" term originating from the cosine of the dual 
field: 

H_(x) = -v b (xr9 x Xr ~ Xl^xXl) 
+™ (XrXl ~ XlXr) + 2im (XrXl ~ XlXr) (17) 
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We introduce two Major ana fields 

The Hamiltonian H_ then describes two massive Majorana fermions, and 
Pr,Li with masses — m and 3m, respectively: 

F_ = F_ m [C]+# 3 m[p] (19) 

Now we observe that = £, £ 2 = rj and £ 3 = £ form a triplet of Majorana 
fields with the same modulus of mass m. There is one more field p with a larger 
mass, 3m. So, 

H = H m [t} + H 3m [p] (20) 

with 

H m [g\ = E {" V ( ® ~ ^ + (-l)^ 3 i™ (21) 

a=l,2,3 / 

The exact degeneracy in the masses and, as a result, 0(3) invariance of the 
Hamiltonian H[£\ is recovered by a duality transformation £ R — ► — £^ (with all 
the other Majorana fields intact) that effectively changes the sign of the mass of 
the third Majorana component. Let us denote the transformed triplet of fields 
byC 

C = (-i) 5a3 (-i)^C- 

The resulting model H m [£\ was suggested as a description of the S = 1 Heisen- 
berg chain by Tsvelik (|14]]). This equivalence follows from the fact that, in the 
continuum limit, the integrable S = 1 chain with the Hamiltonian 

H = Y,[(SnS n+ i) - (Sj n+1 ) 2 ] (22) 

n 

is described by the critical Wess-Zumino model on the SU(2) group at the level 
k = 2, and the latter is in turn equavalent to the model of three massless Majorana 
fermions, as follows from the comparison of conformal charges of the correspond- 
ing theories: 

riwzw _ ^ _ or* 

^ SU(2),k=2 — 2 — ° [ ~ , Major.ferrnion 

The k = 2 level, SU(2) currents expressed in terms of the fields £ a are given by 

Ir,l = ~2 e<lbC £ R , L ^ R > L (^) 
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When small deviations from criticality are considered, no single-ion anisotropy 
(~ D(S Z ) 2 , S = 1) is allowed to appear due to the original SU(2) symmetry of 
the problem. So, the mass term in (^1|) turns out to be the only allowed relevant 
perturbation to the critical SU(2), k = 2 WZW model. 

Thus, the fields £ a describe triplet excitations related to the effective spin- 
1 chain. Remarkably, completely decoupled from them are singlet excitations 
described in terms of the field p. Another feature is that this picture is valid for 
any sign of J±, in agreement with the effective lattice models (|J) and (|1(]) which 
we actually are dealing with. 

Since the spectrum of the system is massive, the role of the so far neglected 
(marginal) part of the interchain coupling (g) is exhausted by renormalization 
of the masses and velocity. Neglecting the latter effect, this interaction can be 
shown to have the following invariant form 

H marg = ^J ± a J dx [ (I%I a L ) + (e R e L ) ( PrPl ) ] (24) 

which, after transforming back from £ a to £ a , reads 

H marg = ±j ±ao J dx[ (e R e L ) - (&i + e R e L ) 

+(e R e L +m-e R e L )(PrPL)) (25) 

In a theory of N massive Majorana fermions, with masses m a (a = 1,2, ...,N) 
and a weak four-fermion interaction 



Hint = \j29a b J dx(C R C L ) , (g ab = g ba ) 

a^=b 

renormalized masses m a estimated in the first order in g are given by 

9ab n A 

'o = "'./ • 2^ 



m « + 7T— m b In i r (26) 

^ 2ttv \m b \ 

Using (|25|) and (^), we find renormalized values of the masses of the triplet and 
singlet excitations: 

M trip = mfl + ^lnA") (27) 
y Arcv \m\ J 

M sing = 3m(l + ^ln^-) (28) 
\ Attv \m\ 
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3 Correlation functions for the identical 
chains. 



Since the singlet excitation with m s = 3m does not carry spin, its operators do 
not contribute to the slow components of the total magnetization. The latter is 



expressed in terms of the k = 2 577(2) currents fl23l) : 

M a ~ I a R + I a L . (29) 

Therefore the two point correlation function of spin densities at small wave vectors 
(|g| << 7r/ao) is given by the simple fermionic loop. A simple calculation gives 
the following expression for its imaginary part: 

where s 2 = uo 2 — v 2 q 2 . Thus the dynamical magnetic susceptibility at small wave 
vectors has a threshold at 2m. 

It turns out that it is possible to calculate exactly the two-point correlation 
functions of the staggered magnetization. This is due to the fact that the corre- 
sponding operators of the Heisenberg chains are related (in the continuum limit) 
to the order and disorder parameter fields of 2d Ising models fl5 |, ||16|| ; the corre- 
lation functions of the latter operators are known exactly even out of criticality 



Using formulas ( |105| ) of the Appendix A, the components of the total (n^ + ^ = 
ni + 112) and relative (n^ = ni — 112) staggered magnetization can be represented 
as 



~ cos \fnO + cos "\A"6L) n ^ ~ sin \fn6 + sin \fnQ_ 

~ sin \A"#+ cos , ~ cos \/tt6 + sin \fn6- 

n z ~ cos V^4>+ cos V^4>- } n l~' ) ~ sin v / tt0+ sin v^r^- (31) 

The fields <f)+,0 + and are governed by the Hamiltonians (0) and @, re- 

spectively. Let us first consider exponentials exp(±i A /7T0 + ), exp(±iy/7c9 + ). Their 
correlation functions have been extensively studied in the context of the noncriti- 
cal Ising model (see, for example, Ref. [fT^|). It has been shown that these bosonic 
exponents with scaling dimension 1/8 are expressed in terms of the order (a) and 
disorder (//) parameters of two Ising models as follows: 

COs(v / VT0+) = <XlO"2, Sin(v / VT0+) = /il/U 2 , 

cos(v / tt^+) = Aiicr 2 , sin(v / 7r^+) = 0^2 (32) 
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Let us briefly comment on this correspondence. 

As already discussed, the j3 2 = An sine-Gordon model H + , Eq.(^), is equiv- 
alent to a model of two degenerate massive Majorana fermions, Eqs.(|H]), (|15|). 
As is well known (see, e.g., Ref. |18|]), a theory of massive Majorana fermion 
describes long-distance properties of 2d Ising model, the fermionic mass being 
proportional to m ~ t = (T — T c )/T c . So, H + is equivalent to two decoupled 2d 
Ising models. Let Oj and fij (j = 1,2) be the corresponding order and disorder 
parameters. At criticality (zero fermionic mass), four products 0"i<7 2 , /ii/i2, cri^2 
and /ii<72 have the same critical dimension 1/8 as that of the bosonic exponentials 
exp(±iA/7T0 + ), exp(±i v /7r^ + ). Therefore there must be some correspondence be- 
tween the two groups of four operators which should also hold at small deviations 
from criticality. To find this correspondence, notice that, as follows from (|7|), at 
m < < cos y/Tt+cj) >t^ 0, while < sin y/Tx<p + >= 0- Since the case m < corre- 
sponds to the ordered phase of the Ising systems (t < 0), < <j\ > = < a 2 > 7^ 0, 
while < Hi > = < ji2 > = 0. At m > (disordered Ising systems, t > 0) the sit- 
uation is inverted: < cos y/7i+(f) >^ 0, < sin a/7t0+ >^ 0, < o\ > = < cr 2 >= 0, 
<yUi> = < / U2>7^0. This explains the first two formulas of Eq.(P2"|). 

Clearly, the exponentials of the dual field 8 + must be expressed in terms of 
(J1/X2 and [i\02- To find the correct correspondence, one has to take into account 
the fact that a local product of the order and disorder operators of a single Ising 
model results in the Majorana fermion operator, i.e. 

~ cos \fn{4> + + 9 + ) ~ 01/ii, £ 2 ~ sin \fn{4> + + 6+) ~ a 2 ^2- 

This leads to the last two formulas of Eq. ( j32[) . 

To derive similar expressions for the exponents of 0_ and 6L, the following 
facts should be taken into account: (i) the Hamiltonian (|17D describing 
modes is diagonalized by the same transformation ([18]) as the Hamiltonian ( |T2"D 
responsible for the "+"-modes; (ii) the Majorana fermions now have different 
masses, and (iii) one fermionic branch has a negative mass. In order to take a 
proper account of these facts one should recall the following. 

(a) A negative mass means that we are below the transition. 

(b) It follows from (ii) that "— " bosonic exponents are also expressed in terms 
of order and disorder parameters of two Ising models, the latter, however, being 
characterized by different £'s. We denote these operators as 03, ^ (mass —m) 
and a, \i (mass 3m). 

(c) Operators corresponding to a negative mass can be rewritten in terms of the 
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ones with the positive mass using the Kramers- Wannier duality transformation 



t — ► —t, a — > n — > a (33) 

Taking these facts into account we get the following expressions for the "— " - 
bosonic exponents: 

COS(V / VT0-) = ^3Cr, Sm(y/lT(f)_) = a 3 /i, COs(v / 7T^) = O^CX, SlD.(y/irO_) = fM 3 fji (34) 



Combining Eqs. (|3~2"D and (|34), from (l3ll) we get the following, manifestly 
SU(2) invariant, expressions: 

u y - ai/i2cr 3 (T, n+ ~ CTicr 2 //30- (35) 

-2, n~ ~ /ii/i 2 o" 3 /i (36) 

It is instructive to compare them with two possible representations for the stag- 



gered magnetization operators for the S = 1 Heisenberg chain ( |14| ): 



S x ~ /iier 2 (73, S' 2 ' ~ <Ti/i 2 cr 3 , S* 2 ~ <Ticr 2/ u 3 (37) 

or 

~ ai/i 2 /i 3 , 5* y ~ /iia 2 yu 3 , S 12 ~ ^1/120-3 (38) 

Agreement is achieved if the singlet excitation band is formally shifted to infinity. 
This implies substitutions a ~ (er) 7^ 0, \i ~ (//) ~ for ferromagnetic interchain 
coupling (m ~ Jj_ < 0), or a ~ (cr) ~ 0, ~ (/i) 7^ for antiferromagnetic 
interchain coupling (m ~ J± > 0). Thus, we observe that, as expected, for 
ferromagnetic (antiferromagnetic) interchain interaction the staggered S=l mag- 
netization is determined by the total (relative) staggered magnetization of the 
two-chain system. 

A more precise meaning of this approximation becomes apparent when one 
considers asymptotic behaviour of the corresponding two-point correlation func- 
tions in the two limits r — > and r — > oo((19|). In the limit r — > 00 they are as 
follows; 

MrK(O)) = G a (f) = ^K (f)+O(e- 3f ) (39) 

<MrK(0)> = G,(r) 
= A, {l + -L [ r ~2 (^(f) - ^ 2 (r~)) - fifo(r)ifi(#0 + \ K l^)} } + °^') ( 40 ) 

where f = rM (M = m or 3m), ^! is a nonuniversal parameter, and it has 
been assumed that M is positive. If M is negative the correlation functions are 
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obtained by simply interchanging a and /i, and putting M —>■ —M (the duality 
transformation Q3"B|)). Therefore, as might be expected, at large distances, a 
difference between the ladder and the S = 1 chains appears only in exp(— 3mr)- 
terms due to the contribution of the excitation branch with M = 3m absent in 
the S = 1 chain. 

In the limit f — > the correlation functions are of power law form: 

G a (r) = G,(f) = # (41) 

T 4 

plus non-singular terms. The ratio of the constants Ai and A 2 is a universal 
quantity involving Glaisher's constant (A): 

M = 2 -i/ 6 A - 3 I (42) 
A 1 P 4 v ; 

A = 1.282427129... (43) 

We conclude this Section by writing down the exact expression for the stag- 
gered magnetization two-point correlation functions. The correlation function for 
spins on the same chain is given by 

(n?(r, xK(0, 0)) = G 2 (mr) (mr)G a (3mr ) + Gj(mr)G (T (mr)G At (3mr) (44) 

The interesting asymptotics are 

(n?(r, xK(0, 0)) = — Z at mr « 1 (45) 
2nr 

in 2 

—ZK (mr){l H — -[(mr) 2 (iS' 1 2 (mr) — if 2 (mr)) — mrK (mr)Ki(mr) 

+^ 2 (mr)]} + 0(e- 5mr ); mr » 1 (46) 

where r 2 = r 2 + f 2 x 2 and 

Z 2 4 / 3 e 



, A sa 0.264 (47) 

The complete expressions for the functions G a ^(f) are given in Ref.([|T9|). For 
the interchain correlation function we get 

K(r, xK(0, 0)) = G\ (mr ) G^ (mr ) ( 3mr ) - G 2 (mr)G (7 (mr)G^(3mr) (48) 

At mr << 1 it decays as (mr)~ 2 ; the leading asymptotics at mr » 1 is the 
same as (|46|) (up to the -1 factor). The difference appears only in terms of order 
of exp(— bmr). The important point is that at mr » 1 the contribution from 
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the singlet excitation appears only in the fifth order in exp(— mr). Therefore it 
is unobservable by neutron scattering at energies below 5m. 

Using the above expressions we can calculate the imaginary part of the dy- 
namical spin susceptibility in two different regimes. For \tt — q\ << 1 we have 

2 cosg_i_[pr<5(c<; — \/v 2 q 2 + m 2 ) + F(u, q)] uj < 5m 

^smx {R '(uj,TC - q;q ± ) = Z{ „ \ n2 fi4 r (49) 

A v ' (l + cosgx) ,°- 264 uj » 5m v 1 

I W Ld 2 — v 2 q 2 

where the transverse "momentum" q± takes values and tt. The factor Z is 
assumed to be m-independent so that at m — > we reproduce the susceptibility 
of non- interacting chains. We have calculated the function F(u,q) only near the 
3m threshold where it is equal to 

144 / 

F(u, q) w r Vcu 2 — v 2 q 2 — 9m 2 (50) 

TTin 2 v 

For |g| << 1 we have 

QmX (B V, T, Q±) = [I + cos 2 (g ± /2)]/(s, m) (51) 
where f(s,m) is given by Eq. (|30|) . 

4 Inequivalent chains; non-Abelian bosoniza- 
tion. 

In this section we consider two interacting spin S = 1/2 chains with different 
intrachain exchange integrals Jij ^ Jj?. It turns out that the most adequate 
approach in this case is non-Abelian bosonization. The reason for this is that 
non-Abelian bosonization explicitly preserves the SU(2) symmetry present in the 
Hamiltonian. The Abelian bosonization approach which does not respect this 
symmetry encounters difficulties. 

As shown by Affleck [2~Cfl, by a mapping from a fermionic theory, the S — \ 



Heisenberg antiferromagnet can be described by a k — 1, SU(2) Wess-Zumino- 
Witten (WZW) model with the following action; 

S k = kW(g) 

w( g ) = -±-Jrr (d.g+d^g) d 2 x + r ( g ) 

T(g) = ^/rf 3 Xe^rr(g + 9 Q gg + ^gg + 5 7 g) (52) 
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where matrix g G SU(2). There can in general be marginally irrelevant pertur- 
bations to this theory, which generate logarithmic corrections to the correlation 
function exponents, but do not change their qualitative behaviour (i.e. power 
law). In general this model describes not just the spin S = 1/2 Heisenberg chain, 
but any (l+l)-dimensional system of fermions with the charge degree of freedom 
frozen out and no gap in the spin sector. 

The WZW model may look unfamiliar, but it is not so difficult to deal with 
since its operators and their exact correlation functions are already known from 
the application of conformal field theory pip (see also p2|). As we have mentioned 



above, the great advantage of the WZW model is that it explicitly possesses the 
SU(2) x SU(2) symmetry of the massless fermion spin sector, and is also critical 
(massless). In (l+l)-dimensions the distinction between relativistic fermions and 
bosons is illusory; one can choose to think about a system in either representation 
(this has been known for some time; hence the "Luttinger liquid"). The WZW 
model is therefore just a way of thinking about the spin sector in terms of bosons; 
just as in Abelian bosonization, one can represent operators from the fermionic 
theory in terms of those of the bosonic theory and vice- versa. From a practical 
point of view, these relations between the two sets of operators can be thought 
of as readymade tools. It is not necessary to worry about their slightly exotic 
appearance or their justification in order to apply them. (But those seeking a 
deeper appreciation are referred to the papers cited above.) 

The bosonized expression for the spin operator of the Heisenberg chain is given 

by0; ^ ^ ^ 

~Su = ^Jr + ^Jl + const(-l) n Tr(g + ^ - g^) (53) 
where the currents are given by; 

J a R = -J-Tr(<9_g)g + T" 

J a L = ^Trg+d+gT" (54) 

(T a are the Pauli matrices - generators of the SU(2) group). These currents 
satisfy the SU(2) Kac-Moody algebra described in the Appendix A. 

Consider two Heisenberg chains coupled by an antiferromagnetic nearest neigh- 
bour interaction. It can be represented like this; 



S = W 1 (h) + W 2 ( 



G r+ G . 



+ X 2 Tr [(g-g + ) ~a] Tr [(h- h+) ~a] (55) 
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where the dynamics of one chain is represented by the matrix g and the currents 
G R t L and the other by h and H r^. The indices 1,2 distinguish between different 
spin wave velocities. Without a loss of generality we can put v\ > v<i- 

The currents have conformal dimensions (1,0) and (0,1); using the formula 
for the conformal dimensions of the matrices for the SU(n) group derived in Ref. 
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n 2 ~l 

A = — -, tt 56 

2n(n + k) v ; 

we get that for n — 2, k = 1, g and h both have conformal dimensions (|, |). 
The A2 term is therefore the relevant interaction, whereas the current couplings 
are only marginal. For this reason, the current interaction will be neglected at 
this stage. Then the interaction can be written as; 

Tr[(g-g+)^]-Tr[(h-h+)^] 
= i {Tr[(g - g+)(h - h+)] - Tr[(g - g + )]Tr[(h - h+)]} (57) 

Making the substitution a = gh + , which leaves the measure invariant, and using 
the remarkable identity f21 |: 



W(ah + ) = W{a) + W(h) + — [ Tra + d-ah + d+hd 2 x (5£ 

2ix J 

we arrive at the following expression for the action: 

S =[W 1 (h) + W 2 (h)] 
— I Tia + d_ah + d + hd 2 x 
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+Wi(a) + A 2 [Tr(a + a + ) - Tr(a + h + + H.c.) + Tr(h + - h)Tr(h + a + - ah)J(59) 

(here d± = |(<9 T =F id x )). 

The identity fl5"8|) is nothing very mysterious. It is simply a generalisation of 
an identity familiar from Abelian bosonization. To see this, consider substituting 
explicitly for the special case of Abelian bosonization, the U(l) fields e 1 ^ 1 and 
e ip02 f or matrices a and h respectively. Then the WZW action, W(a) reduces 
to the action for free scalar bosons as we would expect; 

W (a = e 1 ^) = ^-J d+fad-fadPx =^J [(d x <Pi? + (d T <Pi) 2 } d 2 x (60) 
and the interaction term in (|58[) becomes; 

/ Trr>+<9 nh + d, hd 2 x = — 

2ir 



/ Tra + d_ah + d + hd 2 x = [ (d+fad-fc + d+fcd-fa) d 2 x (61) 
2n J An J 
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The field ah + is e 1 ^ 1 ^ = e l ^~ , and so the identity ([58]) becomes; 

J d + cj)id^(j)id 2 x + j d + (j) 2 d-(j) 2 d 2 x 



Aix 



(62) 



(<9 + </>ifL0 2 + d + <p 2 d-(t)i) <rx 

Therefore the identity (|58|) is just an analogue of the following simple statement: 

(V(0! - 2 )) 2 = (V0i) 2 + (V0 2 ) 2 - 2V0i • V0 2 (63) 

where the last term is the "interaction term" . 

We shall consider the most relevant interaction, Tr(a + a + ) first. The effective 
action for a is in this approximation; 

S = W^a) + XTt (a + a + ) (64) 

From the first order RG equation we get 

dX , 1 . , , . 

(2--)A (65) 



dlnL y 2' 

Integrating up to a scale where the coupling becomes of order 1 and taking this to 

2 

give some estimate of the dynamically generated mass, one gets M ~ As. Much 
more information can be found by realising that the model (Bl) is equivalent to 
the f3 2 = 2tt sine-Gordon model (see for example |20f ). 

Thus on the scale |x| >> M -1 the fluctuations of the a-field are frozen and 
we can approximate 

Tr(ah + )Tr(h) « (Tra) : [Tr(h)] 2 : (66) 

At this large scale the cross term containing derivatives of h and a gives the 
irrelevant contribution 

S int ~ M~ 2 d + d^h + d + d_h (67) 

Therefore the asymptotic behaviour at large distances is governed by the fol- 
lowing action: 

S = W x {h) + W 2 (h) + c 2 : [Tr(h)] 2 : (68) 
where c 2 ~ A 4 / 3 and which can be further modified by the coordinate rescaling: 

Xq = y/ViV 2 T, Xi= X (69) 
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such that we finally have 



S = S + S 1 (70) 
So = ^rj d 2 xTr (<9 M h + <9 M h) d 2 x + 2r (h) (71) 

Sx= f d 2 xc 2 {: [Trh 2 ] : + : [Tr(h+) 2 ] : - : [Tr(h - h+)] 2 :} (72) 

where 

1/c = \Jvi/v 2 + \Jv 2 /vi 

The model with action ([70]) is not critical; coupling constants Ci,c 2 undergo 
further renormalization. Let us show that the coupling c 2 renormalizes faster 
to strong coupling. To show this we shall suppose that this is the case and check 
that the obtained result is self-consistent. It is easy to check that the effective 
potential (|72|) vanishes if h is a traceless matrix and has a fixed determinant: 

h « i(an), n 2 = 1 (73) 

Excitations, which correspond to configurations where Trh ^ 0, acquire a gap. 
The estimate for this gap is 

Ml ~ cy/v^ c 2 ~ min(f 1; v 2 )X 4/3 ~ ^M 2 (74) 

On energies smaller than the gap one can treat the h-matrix as traceless. Sub- 
stituting expression (|73|) into Eq. (|70|) we get the 0(3)-nonlinear sigma model as 
an effective action for small energies: 

S = ^Jd 2 x(d fl n) 2 , n 2 = l (75) 

I = {yfcfa + yfefa)(l - (n 2 )) (76) 

The reason why the Wess-Zumino term effectively disappears from the action is 
the following. After substituting Eq.(^) into the expression for T(h) the Wess- 
Zumino term reduces to the topological term: 

2T{ian) = l - J d 2 xe^ u {n^n x d v n\) = 2mk (77) 

where k is an integer number. The factor in front of the topological term is such 
that its contribution to the action is always a factor of 27ri and therefore does not 
affect the partition function. The mass gap of the model ( [75] ) is given by 

M = Mod' 1 exp(-27r/5) w M[l - ^- \n{M/M )\ exp(-27r/c) (78) 
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As long as this gap is much smaller than M , the adopted approach is self- 
consistent. The latter is achieved for any appreciable difference between the 
velocities. 

Excitations of the 0(3)-nonlinear sigma model are S = 1 triplets (p3"l). Thus, 
the spectrum is qualitatively the same as for identical chains. That is what one 
might expect because the model of Majorana fermions is a strong coupling limit 
of the 0(3)-nonlinear sigma model (see Ref. |T4|j). 

The correlation functions of the 0(3)-nonlinear sigma model are known only 
in the form of the Lehmann expansion (|[24||). 



(n(r, x)n(0, 0)) ~ K {mr) + 0(exp(-3mr)) (79) 

Note that the first term in the expansion coincides with the one for identical 
chains. Therefore a difference in dynamical magnetic susceptibilities for both 
cases will become manifest only at energies u > 3m. The lowest feature in 
9mx' R ' (uj , q) is in both cases the sharp peak 

m 

^Jq 2 + m 2 



9mx (fl) (oj, q) ~ 6{u - Jq 2 + m 2 ) (80) 



corresponding to the triplet excitation. Such a peak has been observed in 

(vo) 2 p 2 o 7 ii|. 



5 String order parameter in the spin-ladder 
model 

Den Nijs and Rommelse p5f (see also p6| ) have argued that the gapful Haldane 
phase of the S = 1 spin chain is characterized by a topological order measured 
by the string order parameter 

m— 1 

<O a >= lim < S>xp(i7T S D S ™ >. (^=1. a = x,y,z) (81) 

1 1 j=n+l 

The nonzero value of < O a > has been related to the breakdown of a hidden 
Z 2 x Z 2 symmetry P7| . In this section we use the Abelian bosonization method 
(section 2) to construct the string operator in the continuum limit of the S = 1/2 
spin-ladder model and identify the corresponding discrete symmetry with that of 
the related Ising models. 

Since spin-rotational invariance remains unbroken, the string order parameter 
must respect this symmetry. However, Abelian bosonization is not an explicitly 
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SU(2) invariant procedure. For this reason, it turns out that it is the ^-component 
of the string operator that acquires a simple form in the continuum limit. On the 
other hand, due to the unbroken SU(2) symmetry, the very choice of the quanti- 
zation (z-) axis is arbitrary; therefore the expectation values for all components 
of the string operator will coincide. 

To construct a string order parameter O z (n,m) for the spin-ladder model, we 
shall follow the same route as that previously used for the bond-alternating S = 
1/2 chain |27j (technical details are given in the Appendix B.). We start from 
the lattice version of the model, construct a product of two spin-1/2 operators 
belonging to the j-th rung, S( (j^fC/), and then take a product over all rungs 
between j = n and j — m: 

m I m \ 

O z (n,m) = U(-4SZU)SZ(j)) = ex P i7r EI S tU) + %U) ] (82) 

j=n \ j=n J 

Assuming that \m — n\ ^> 1, we pass to the continuum limit in the exponential 
and retain only the smooth parts of the spin operators expressing them in terms 
of the spin currents J z . RL (x), (a = 1, 2): 



O z (x,y) 




) + JUlV)] ■ (83) 



Using Eqs.fl5D,(|97D, we find that the exponential is expressed in terms of the field 
cf) + only. Thus we find a very transparent representation for the string operator: 

O g (x, y) = exp(iv^[0 + (x) - 0+(t/)]) (84) 

Using Eq.(0), 

exp (i^(J) + (x)) ~ ai<j 2 + i^i^2 (85) 

we find that the string operator is expressed in terms of the Ising order and 
disorder operators. For either sign of Jj_, we find that, in the limit \x — x'\ — > oo, 
the vacuum expectation value of O z (x,y) is indeed nonzero: 

lim (O z (x,y)) ~ < a x > 2 < a 2 > 2 = < a >V 0, J± < (86) 

\x— x'\— >oo 

lim (O z {x,y)) ~ < ^ > 2 < /x 2 > 2 = < fi >V 0, J ± > (87) 

\x— x'\— >oo 
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As in the case of the bond-alternating spin chain, the nonvanishing expectation 
value of the string order parameter in the limit of infinite string manifests break- 
down of a discrete Z 2 x Z 2 symmetry. This is the symmetry of two decoupled Ising 
models described by the Hamiltonian H + in the Majorana fermion representation 
(|14|): H + = Hml^ 1 ] + H m [^ 2 } remains invariant with respect to sign inversion of 
both chiral components of each Majorana spinor, £% L — > — £% L , (a = 1,2.) Un- 
der these transformations, the Ising order and disorder parameters change their 
signs. On the other hand, since the two Majorana fermions are massive, this 
symmetry is broken in the ground state of H + : the mass terms break the duality 
symmetry — > — ££, £^ — > £r- This amounts to finite expectation values of the 
Ising variables G\ and 02 (or \i\ and /Z2), which in turn results in a nonzero string 
order parameter, as shown in Eqs.(p6[) and (|87|). 

6 Conclusions 

As the reader might see the spin ladder presents an exciting opportunity to study 
a formation of massive spin S = 1 and S = particles which appear as bound 
states of spin S = 1/2 excitations of individual Heisenberg chains. At small 
interchain coupling | J±\ << J\\ the masses of these particles are of order of \J±\. 
The S = 1 branch is always lower independently of the sign of J±. At J±/ J\\ — > 
the singlet spectral gap is three times as large as the triplet one. The comparative 
smallness of the interchain coupling creates room for the crossover from the single 
chain behaviour to the strong coupling regime which resembles the S = 1 chain. 
The imaginary part of the dynamical spin susceptibility ~x!'{oj, q; q±) calculated in 
Section 3 contains essential information about this crossover. At small energies 
the susceptibility exhibits a sharp peak around q = ir correspondinding to the 
stable S = 1 massive particle; at energies u > 3m Q) exhibits an incoherent 

tail originating from the multi-particle processes. Below the 5m-threshould the 
singlet branch does not contribute to ^'(uj,g) and the latter coincides with the 
susceptibility of a S = 1 chain. The contribution from the singlet mode becomes 
essential at energies much greater than the spectral gap and the susceptibility 
asymptotically approaches its value for a spin-1/2 chain. We emphasise that the 
described picture holds only in the ideal limit J±/J\\ — > 0. We suppose that in 
real systems it will be difficult to make this ratio less than 0.1. 
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Appendix A: Basic facts about bosonization 



Antiferromagnetic spin chain Hamiltonians, such as the Heisenberg Hamiltonian 

N 

H = j£S n -S n+1 , (5 = 1/2, J>0) (88) 

n=l 

can be mapped onto fermionic theories. Using bosonization, these can be recast 
as generalised Sine-Gordon or WZW models. This is useful because a great deal 
is known about these theories, such as correlation functions, scaling dimensions 
of operators etc. A brief summary of this approach is given below. 

Following Refs. ||20[| , we start from a symmetry preserving fermionization of 
the spin operators 

Sn = ^rJf^nP (89) 

To eliminate the redundant zero- and double-occupancy states, the constraint 
J2 a i'na'tpna = 1 for all lattice sites n should be imposed. Such a constraint will 
effectively work, if one considers a 1/2-filled, U > Hubbard model for the field 
ipna- In this model, a Mott-Hubbard charge gap m c is known to exist for any 
positive U. Therefore, at low energies, \E\ <C m c , only spin excitations remain; 
those describe universal dynamical properties of the spin-chain model (|8~8"D in the 
continuum limit. 

Assuming that U <C t, we linearize the free-particle spectrum near two Fermi 
points, ±fcp (kp = 7r/2ao), and decompose the Fermi field into right-moving and 
left- moving chiral components: 

i>na -> V^O 4(4 tpa{x) = (-i)* ^ 'Ra(x) + ^^La(x) (90) 

We then introduce the scalar [U(l)] and vector [SU(2)] currents (the local charge 
and spin densities) 

JR,L = ■ ^R,L;a^R,L;P : ; J R,L = ■ ^ija^V' R,L;f3 ■ (91) 

satisfying anomalous (U(l) and SU(2)) Kac-Moody algebras: 

[Jr(x),J r (x')} = -6'(x-x'), (92) 

17T 

[J a R (x), J b R (x')) = ie abc J R (x)6(x - x') - ^S ab 5'(x - x') (93) 

(with similar relations for the left components). These algebras lead to fermion- 
boson duality which allows to represent the Hamiltonian of free fermions as a 
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sum of two independent (commuting) contributions of gapless charge and spin 
collective modes (Sugawara form): 



n — n V{\) "T U SU(2) 



J dx[ (: JrJr : + : J L J L :) + (: J R ■ J R : + : J L • J L :)] (94) 



^(:J R J R : + :J L J L : ) + *f 

The charge part is equivalently described in terms of a massless scalar field <p c . 
Under identifications Jr + Jl = -^d x (p c , Jr — Jl = — -^n c , where Il c is the 
momentum conjugate to the field C , one obtains 

H° m = jjdx [Tll(x) + {d x ct>c{x)f] (95) 

The spin part Hg V ^ represents the level k = 1 5'L r (2)-symmetric critical Wess- 
Zumino-Witten (WZW) model. 

A weak Hubbard interaction preserves the important property of charge-spin 
separation, Huubbard — H c + H s . Umklapp processes relevant at 1/2-filling trans- 
form to a quantum sine-Gordon model 

H c = Jdx [| (n c 2 + {d x <f) c ) 2 ) + canst g cos &0 C ] (96) 

which at g ~ U/t > occurs in its strong-coupling, massive phase {(3 2 < 8n), with 
the single-soliton mass m c being just the the Mott-Hubbard commensurability 
gap. 

In the spin sector, interaction —2gjR • Jl is added to Hg V , 2 y This interaction 
is marginally irrelevant (since g > 0). Therefore, the universal scaling properties 
of the Heisenberg S = 1/2 spin chain flggp are described by the level k — 1 WZW 
model H° sm 0. 

The possibility of an Abelian bosonization of the Heisenberg chain (|88|) stems 
from the fact that conformal charges of the k = 1 SU(2) WZW models and 
free massless Bose field coincide: C^f(^ k=1 = Cb OS on = 1- Using relations 
(1/3) Jr(l) ■ Jr(l) = Jr(l)Jr(l)i Hsu(2) can b e expressed in terms of J 2 -currents 
only; introducing then a pair of canonical variables, <p s and U s , via 

JR + Jl = ^=d^ si j* _ j- = __ ^n., (97) 

V27T V27T 



one finds 

H° su{2) ^H B = V ^Jdx [U 2 s (x) + {d x cj> s {x)f) (9f 
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The price we pay for this simplification is the loss of spin rotational invari- 
ance in the bosonized structure of the spin currents: the J x and J y cannot be 
represented as simply as J z , and require bosonization of the Fermi fields: 

^R,L-Ax) ~ (2na )- 1/2 exp (±i^ (p R>L ;a(x)) (99) 

Linear combinations 

constitute scalar fields $ a and their dual counterparts Q a introduced for each 
spin component. The fields describing the charge and spin degrees of freedom are 
defined as follows: 



where d x 9 CtS = n cs . 

T± 



(100) 



To bosonize J RL , use fl99|) to obtain: 



j r = ^R]^Ri = exp(-zv / 27r(0 s - 9 S )) 

Jt = ^Li = -^—exp(iV^((l> a + 9 a )) (101) 

Note that, as expected, the charge field C does not contribute to the spin SU(2) 
currents. Moreover, despite the fact that the definitions ( |101|) contain cutoff ao 
explicitly, the current- current correlation functions are cutoff independent and 
reveal the underlying SU(2) symmetry: 

s:ab i 

< .r{x)j\x') >= - — - (102) 

Aty z (x — x y 

The SU(2) currents J r(x) , J l(x) determine the smooth parts of the spin op- 
erators in the continuum limit. Namely, at ao — > 

S n -> a S(x), S(x) = 3 R (x) + J L (x) + (-l) n n(x) (103) 

where 

"(*) = A(*)y^^)+^ (104) 
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is the staggered part of the local spin density. 

When bosonizing ( |104 ), the (redundant) charge excitations emerge, since off- 



diagonal bilinears like ip^ipL and ip^ipR describe particle-hole charge excitations 
with momentum transfer ±2kp- We find: 

n z = sin y/27i(f) c cos y/27i(f) s 

Too 

= sin(v / 27T0 c /2) exp(±iy/2n9 s ) 

Being interested in the energy range \E\ <C m c , one can replace the charge op- 
erator sin(\/27r</> c /2) by its nonzero vacuum expectation value; we denote this 
(nonuniversal) value by A = — < sin(v^27T0 c /2) > and arrive at bosonization 
formulas for n(x): 

n z {x) = cosy/2n(j) s (x) 

■ri^tx) = — exp[±iV2^9 s (x)] (105) 

This completes the bosonization of the spin operators for the isotropic Heisen- 
berg chain. Notice that the critical dimensions of the smooth and staggered parts 
of the spin densities are different: 

dimJ a = l, dim n a = 1/2. (106) 



Appendix B. Hidden Z2 x Z2 symmetry and 
string order parameter in the bond-alternating 
S = 1/2 Heisenberg chain 

In addition to the S = 1/2 spin-ladder model, there is another system which is 
related to the S = 1 spin chain - the spin- 1/2 chain with alternating ferromagnetic 
and antiferromagnetic bonds: 

N/2 

H = 4/£[(S 2i _ 1 ■ S 2i ) - /?(S 2j • S 2i+1 ] (107) 

This model is instructive in the sense that the string order parameter, whose 
nonzero expectation value signals breakdown of a hidden discrete symmetry, can 
be easily constructed [27]. The analogous construction is then directly generalized 



for the spin-ladder model. 
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A gap in the excitation spectrum of the model (|107|) persists in the whole range 
< [3 < oo. At (3 — the ground state of model represents an array of discon- 
nected singlets. At (3 ^> 1, strong ferromagnetic coupling between the spins on the 
< 2j, 2j + 1 >-bonds leads to the formation of local triplets, and the model (|107| ) 
reduces to a S = 1 Heisenberg chain. Using a nonlocal unitary transformation, 
Kohmoto and Tasakif^] have demonstrated equivalence of the model (|107|) to a 
system of two coupled quantum Ising chains, i. e. two coupled 2d Ising models. 
This transformation provides an exact representation of the spin operators as 
products of two Ising-like order (cr, r) and disorder (a, f ) operators, essentially a 
lattice version of relations (p2[), (|34l) (see, e.g., p8| ). Nearest-neighbor bilinears 
of the original spin operators take the form 



where 



A QX QX Z Z 

w 2j°2j+l — °j°j+l> 


4^- 


QX 

-1^2j - - 


T X 
3 


aqV qy _ z z 
^°2j°2j+l — T j T j+ll 




qy 

- 


a 3 


A QZ qZ Z Z 2 2 

^j^j+l ~ a j a j+l T j T j+l' 


4^- 


qz _ _ 


°3 T 3 






1/2 T j+l/2 




N/2 




3-1 





(106 



(109) 



^ 2 + i/2= n °f> ^ + i/2 = n^ (no) 

1=1 



Relations (108) make the Hamiltonian (|107|) equivalent to two coupled quan- 



tum Ising chains: 

N/2 
3=1 



The model ( 11 1|) is invariant under independent rotations of the a and r spins by 



angle ir about the spin x-axis which comprise a Z 2 x Z 2 group. Since this group 
is discrete, it can be spontaneously broken, in which case the spectrum of the 
system would be massive. It is easily understood from (|111 ) that, in the limit of 



large positive (3 when the model reduces to the S — 1 chain, the Z 2 x Z 2 -symmetry 
is broken, with 

< o] > = < t* > = < o)t* > y£ 0, (112) 
(It has been used in Eq. (|112 ) that, under transformation //| = cr|rj to a new pair 



of variables, pfc and r?, the two-chain Hamiltonian (|111|) preserves its form.) 
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Representation (\108\) hints to the way how an order parameter measuring 
breakdown of the Z 2 x Z 2 -symmetry should be constructed out of the spin oper- 
ators S". Following Kohmoto and Tasaki, consider a product 

2n— 1 ri— 1 n—1 

n asf = n = nc-^+i) 

l=2k j=k j=k 

= (-i) n ^(^^ +1 )K +1 ^ +2 ) • • • = (-irvx (us) 

Using the relation icr" = exp(27rcr"/2), we find that 

/ 2n-l \ 

<?*(*, n) = exp Utt £ 5f = (H4) 

\ l=2k J 

This is the x-component of the string order operator. According to ( |112| ), in the 
limit | A; — n\ — > oo, its vacuum expectation value is nonzero: 

< O x {k,n) > -> < a >V (115) 

It is important that the string always contains an even number of sites, starting 
at an even site and ending at an odd site. For a string starting at an odd site and 
ending at an even site, the corresponding string operator is expressed in terms of 
disorder operators and therefore has zero expectation value: 

2n n 

II 2Sf = Y\_ ^2j-1^2j = ( — 1)™ (^fc+l/2^fc+3/2)(^fc+3/2^fc+5/2) ' ' ' (^n-l/2^n+l/2) 



/=2fe+l j=k+l 



i \n~k~z ~z 

L ) ' k+1/2 1 n+1/2 



The y and z components of the string operator are constructed in a similar 
manner: 

(2n-l \ 2n-l n-1 
i7r e s ?) = n 2i5? = nMsiAVi) 
/=2fc / «=2fc j=k 

n-1 

= n = ^ ( iie ) 

j=fe 

(2n-l \ 
iw E 5? }=°kTk« (H7) 
Z=2fc / 

The SU(2) invariance of the expectation value of the string order parameter 

(2n-l \ 
ivr ^ 5f] , (5 = 1/2, a = x,y,z) (118) 
Z=2fc / 

follows from ( |112| ). 

Notice that in the limiting case /3 > 1, the string order parameter ( |118| ) for 
the S = 1/2 bond-alternating chain automatically transforms to the exponential 
of the string order parameter (jSll) for the S = 1 chain. 
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